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Abstract: We consider super symmetric QCD in the free magnetic phase with massless 
and massive flavors. The theory has a supersymmetry breaking pseudo- moduli space of 
vacua and a runaway behavior far away from the origin. A two-loop computation reveals 
that the origin is destabilized and there is no meta-stable SUSY breaking solution. We also 
study the embedding of this model in type IIA string theory and find evidence for similar 
behavior. The perturbative brane dynamics involves simple interactions between branes, 
correctly predicting the two-loop result in the gauge theory. Our results also apply to the 
case when all the flavors are massive but have hierarchy among them, leading to possible 
instability which is manifest both in field theory and the brane description. 
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1. Introduction 

Dynamical SUSY breaking (DSB) may be an attractive explanation of the hierarchy be- 
tween the electro- weak scale and the Planck scale Q. In spite of the fact that there are 
models which break SUSY spontaneously at the ground state (see for example the re- 
view Q), these models are extremely rare and they are subject to some severe constrains. 
Moreover, a generic (and calculable) SUSY breaking model should have a spontaneously 
broken R-symmetry ||, leading to a relatively light R-axion, which gets its mass only from 
couplings to supergravity § or to some higher dimensional operators. Current astrophys- 
ical observations disfavor this possibility. 

These difficulties provide a good motivation to consider the possibility that SUSY is 
broken dynamically in a meta-stable vacuum (see Q for a review of DSB both in stable and 
meta-stable vacua). The idea of DSB in a meta-stable vacuum got a lot of attention after 
it was shown by Intriligator, Seiberg and Shih (ISS) in || that a simple and generic class 
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of models, supersymmetric QCD (SQCD) in the free magnetic range, possesses local meta- 
stable SUSY breaking vacua (see e.g. 0, |, |, 0). The authors of § considered SQCD 
with gauge group SU(N C ) and Nf flavors in the range N c < Nf < 3N c /2. It was shown 
that if the quarks are given small and equal masses then the theory has a long-lived SUSY 
breaking vacuum. Since in this range the magnetic description of the theory is weakly 
coupled in the IR, the analysis at low energies was done using the Seiberg duality |TT| . 
At tree-level massive SQCD has a moduli space of SUSY breaking vacua, but one-loop 
quantum effects stabilize the pseudo-moduli at the origin of field space. 

In this work we study SQCD with massless and massive (but light) flavors. There are 
several reasons to consider this model. First, such a theory is an extreme case of massive 
SQCD with generically distributed masses. Even though the ISS model looks as a good 
starting point for direct gauge mediation (for early models see H, ||, [lj, IH)> ii; should 
be modified in order to produce viable phenomenology. One of the reasons for such a 
modification is an approximate R-symmetry at low energies which protects gauginos from 
getting masses. Some of the modifications included hierarchical quark masses |8| since this 
introduces another scale with which one can tune the lifetime independently of gaugino 
masses. Therefore, it is interesting to understand a limiting case of a generic distribution, 
taking some of the flavors to be completely massless. 

Another motivation is that SQCD has a natural embedding in type IIA string theory, 
as the low energy theory on intersecting Neveu-Schwarz (NS) fivebranes and D-branes |l6[ 
(for a review, see fj"7| ). The brane description of the meta-stable vacua of || was studied 



recently in [18, |T^, ^TJ. Many modifications and variations of the basic model of Q 
were constructed, along with their brane descriptions (see e.g. |^, [24]]). It is tantalizing 
that in these examples one could identify identical patterns of meta-stable SUSY breaking 
states in the gauge theory and the classical brane system. In gauge theory, pseudo-moduli 
are stabilized by one-loop quantum effects (f| 0, ^5| , while in the classical brane dynam- 
ics regime, gravitational attraction in the NS fivebrane background stabilizes the branes 



in long-lived SUSY breaking meta-stable configurations [21, 24]. It is interesting to see 
whether this correspondence can be pushed further, and to check whether also our system 
has similar qualitative properties in the perturbative brane dynamics regime. 

We consider SU(N C ) SQCD with NfQ massless flavors and Nf — NfQ massive ones in 
the range < NfQ < N c < Nf < 3N c /2 and study it in the dual magnetic description. In 
such a case the maximal possible rank of the quarks mass matrix is still larger than the 
rank of the dual gauge group, hence, there is no classical supersymmetric solution. Instead, 
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classically, these models possess a moduli space of SUSY breaking vacua. However, the 
pseudo-moduli associated with the massless (electric) quarks are not lifted by one-loop 
quantum effects in field theory p^] and a two-loop calculation is required to decide what 
is the fate of this system. 

In this paper we perform the calculation of the two-loop effective potential for these 
pseudo-moduli. We show that at the two-loop level these directions are destabilized and, 
consequently, there is no SUSY breaking meta-stable vacuum near the origin. This result 
is also important for the case when all the flavors are massive, but there is mass hierarchy 
among them. In that case the two-loop contribution of heavy quarks will dominate the one- 
loop contribution of light quarks and the SUSY breaking solution of may be destabilized. 

In addition, we study the corresponding brane description and find compelling evi- 
dence that a similar instability occurs there. In particular, in the appropriate sense, the 
"origin" is destabilized by the brane dynamics. Note that so far in all the studied exam- 
ples it was found that there is a non-trivial correspondence between the weakly coupled 
brane dynamics and field theory: whenever there is a meta-stable state in gauge theory one 
could identify a meta-stable state in the classical branes picture. Our work provides an- 
other non-trivial check of this correspondence, beyond one-loop effects in field theory (and 
beyond classical gravity in the brane dynamics). We emphasize that understanding the 
perturbative brane dynamics involves simple classical considerations, correctly predicting 
the result of an intricate two-loop evaluation in field theory. 

This paper is organized as follows. In section § we study a simple Wess-Zumino 
(WZ) model which has a similar structure to the low energy SQCD, and then turn on 
the appropriate gauge interactions. In section || we present (after a brief review) the 
brane description of this gauge theory and analyze it. In section ||] we comment on some 
implications of our results to a theory with a general distribution of masses, especially 
the issue of stability of the local minimum of massive SQCD and its brane construction. 
Finally, we summarize in section ||. Appendix |A] contains a brief review of [26| and some 
technical details related to our calculation. 



2. Field Theory Analysis 

2.1 A Simplified Wess-Zumino Model 

In this section we analyze a simple WZ model, which has a pseudo-moduli space of SUSY 
breaking vacua at the one-loop approximation, and show that it does not have any SUSY 
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Table 1: The chiral superfields and their global U(l) charges. 



breaking minimum near the origin of field space. Consider a model with the chiral super- 
fields listed in Table |l[ and a superpotential 



w = V^-- V s ($n + $22), 



(2.1) 



where i, j = 1, 2, 3. The components of the matrix $ and vectors q, q are given in terms of 
the fields in Table || by 



/ $11 $12 x\ 

$21 $22 Y 
\X Y Zj 



X P cr , 



/x\ 

P 



(2.2) 



Let the Kahler potential be canonical. The parameter h controls our loop expansion. Of 
course, the physical parameter corresponding to h is IR free, allowing a faithful perturbative 
treatment. 

The model has manifest SU(2) x U(l) xp x U(l) a x U (1)5- x U(1)r symmetry, where the 
two upper components of q and q transform as fundamentals of the SU(2) and $ transforms 
in the adjoint. 1 Under the various U(l) symmetries the fields transform as summarized in 
Table |. Note that once n is turned off there is an S"f7(3) 2 x U(\)b x U(l) R symmetry. 
The baryon number is still present in our model as U(l) xp + U(l)a — U(l) a = U(1)b- 

This system has no classical SUSY preserving vacuum. The F-terms of the relevant 



1 By that we mean that the upper-left 2x2 submatrix of $ sits in the adjoint, (X,Y) and (X,Y) are 
fundamentals of 517(2) and Z is neutral. 
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meson components are 

dW 



X P 



hfJ,% x2 , i, 3 = 1,2. (2.3) 



* = 0, q = (ji o) , q 



The first term is at most of rank one while the second term is of rank two. Thus, SUSY is 
broken by a rank condition. Nonetheless, there is a stationary point with positive energy 
(i.e. spontaneously broken SUSY), 

/A 

• (2.4) 

w 

The global symmetry is broken as SU(2) x U(l) xp U(l)'. Thus, the above classical 
solution has three Goldstone bosons. In addition, classically, there are some dangerous 
pseudo-flat directions. Our purpose is to understand their quantum mechanical fate. For 
convenience, we take /i and h to be real and define 

P± = -J=(P±P), X± = -/|(X±X)- (2.5) 

The squared mass of the fields a, a, X, X, $12, $21 is h 2 p? . Similarly, the mass of x+i $n> 
3f?p_ is 2h 2 fj, 2 . The three real Goldstone bosons are Qx-^P+^P-- All the rest are 
"accidental" pseudo-moduli which are not protected quantum mechanically, in general. 

The results of the one-loop effective potential in this model are known from [2^] and 
we shall review them here. All pseudo-moduli fields but Z obtain similar positive mass 
squared terms, 

2 2 2 t 4 2 m ^ — 1 

™*R X _ = ™<K 22 = 2m Y ,Y = k V 8?r 2 ■ 

However, Z remains massless at one- loop. One can argue that in the one-loop effective 
potential there will be no Z n terms, for any n > 1. The way to see it is to turn off the 
expectation value of the classical pseudo-moduli Y and Y. Doing so, particles whose mass 
depends on the expectation value of Z are decoupled (in the mass matrix) from the p 
sector which breaks SUSY. Thus, they sit in supersymmetric multiplets and the one-loop 
contribution vanishes identically. 

This means that in order to understand the dynamics of this model it is necessary 
to compute the two-loop effective potential along the pseudo-moduli space parameterized 
by Z. Explicitly, we replace all the fields by their fluctuations and assume, without loss 
of generality, that the field Z obtains a real expectation value around which it fluctuates. 
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Indeed, we can use the symmetry generator U(l) a + U(l)& to rotate the point where Z is 
real to any other complex value of Z with the same magnitude (Note that all the other 
expectation values vanish since they correspond to fields which are massive at tree-level or 
one-loop.). 

The superpotential is given by 



-hp 2 S^ 2 2 - h[i 2 5$n + 



(2.6) 
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+ 2 
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The spectrum of masses is as quoted above (when expanding around Z = 0), except that 
the fields a, a, X, X mix in a simple manner. The mass eigenstates are given by some linear 
combinations 



8a = sin6 5 A + cosfl SB, 



SX = cos9 SA- sinO SB, 



(2.7) 



and analogous equations for the tilded fields (with the same mixing angles). Hereafter 
we use the notations sg = sine and eg = cos 9. A and B are mass eigenstates with the 
following masses 



^ 2 + ^f^Tv 



The mixing angle is 



4 



m 



m 



B 



(2.8) 
(2.9) 



The mass spectrum of all the particles except p± is supersymmetric. 

From now on, the two-loop evaluation is, in principle, straightforward (but in practice 
there are many diagrams). All the required two-loop functions and diagrams are beautifully 



described in [26]; some highlights are reviewed in Appendix A. We have simplified the 
computational task (in particular, the number of diagrams) with a few tricks which may 
be useful also in other models. 

Consider a different theory in which we switch off the linear term for $22 in the 
superpotential W. In other words, we consider a theory whose superpotential is W = 
yV+hfi 2 8&22, where W is given by d2.6|) . In this model the moduli space parameterized by Z 
still exists but now it is a supersymmetric moduli space. It cannot be lifted by perturbative 
quantum corrections. This means that the effective two-loop potential vanishes identically 
as a function of Z. Thus, we can write the trivial equation, 



(2) 



(2) 



(2) 



(2.10) 
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(2) 

for the two-loop effective potential we are after, Vyy . Note that all of the Yukawa, cubic and 
quartic interactions are identical in the two models. In fact, the only difference is that the 
fields p± of the model W' are in supersymmetric multiplets with mass h 2 p 2 . Consequently, 
diagrams that do not cancel on the right hand side of Q2.10 ) contain necessarily a p± scalar. 



However, this is not the only simplification we can make. Since we want the diagrams to 
have some Z dependence, we should better have an A or B (fermion or boson) running in 
the loop. Otherwise, the diagram contributes only to the overall zero-point energy which 
we are not interested in. In this way we remain with only three different diagrams! They 
are depicted in Fig. [I]. 



Y,Y 



/ 



/ 

AS 



A,B 



Figure 1: The only 3 two- loop diagrams contributing to the Z dependent part of the effective 
potential. Conforming with pfJ, we refer to them as SS, SSS and FFS, respectively. 



Note that in general there is another possible topology for a two-loop diagram, one 
that includes mass flips for fermions (the third diagram depicted in Fig. ||). Even though 
the fermionic mass term in our theory is not diagonal this diagram is absent because there 
is always a ipy fermion in the loop, which is massless. 

At this stage, it remains to evaluate the coefficients of the diagrams in Fig. |l[ Of 
course, as follows from fl2TQ), we must subtract from each of the diagrams the correspond- 
ing diagram in the theory W'. In terms of the functions given in and reviewed in 
Appendix A, we get 

y(2) _ y(2) , y(2) y(2) , * 

v — V SS + V SSS + FFS' \ Z - Ll ) 

where 

Vss = h2 4 {fss{2h 2 p\m\) - 2f ss (h 2 p 2 ,m 2 A )) + (A, s 2 ) ~ (B, c 2 ), (2.12) 

V SSS = + Zs e ) 2 (f SS s(0,0,m 2 A ) + fsss{0,2h 2 f/,m 2 A ) - 2f SSS (0, h 2 p 2 , m 2 A )) + 

+ {A,c e ,s e ) <-> (B,-S0,co), 

(2.13) 
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V F % = h 2 s 2 e (f FFS (0,m 2 A ,0) +f FFS (0,m 2 A ,2h 2 n 2 ) -2f FFS (0,m 2 A ,h 2 iJ?)) + 

+ (A,s 2 )~(B,c 2 ). (2.14) 

Our results manifestly look like the difference of amplitudes in two different theories. 
From here on, it is straightforward to expand these functions in a Taylor series and to 
get that the overall contribution to the effective potential is (promoting Z to be complex 
again) 

V® = const + h 6 fi 2 f-l-?L + In 4^ \Z\ 2 + 0(\Z\ 4 ). (2.15) 

Thus, the origin is destabilized. An examination of the effective potential as a function 
of Z shows that there is no minimum around the origin; the effective potential decreases 
monotonically. 

Not surprisingly, there is no dependence on the renormalization scale, Q, in front of 
\Z\ 2 . It is a consequence of the following RGE argument. The effective potential satisfies 
an equation of the schematic form 

( d d d\ 

where (5h is the beta-function of a (physical) coupling h and 7^ is the anomalous dimension 
of a field 4>. Since both of these functions begin at one-loop order (or higher), and since 
there is no \Z\ 2 term at one-loop or at tree-level, Q-g^j should annihilate the two-loop 
coefficient of \Z\ 2 (or higher powers of Z), as indeed happens. Such an RGE argument is 
more general: loosely speaking, this means that whenever a physical effect appears for the 
first time in the effective potential, it must be renormalization scheme independent. 

2.2 Supersymmetric QCD 

Our general model is SQCD, whose UV electric description is given by the superpotential 

N f 

W = zZ m (a)QaQ a , (2.16) 
a=l 

where Q a (Q a ) is in the (anti-)fundamental representation of the gauge group SU(N C ). 
We choose N c < Nf < 3N c /2, where the theory is in the free magnetic phase, and we 
take NfQ of the flavors to be massless, such that < NfQ < N c . This implies that, non- 
perturbatively, far away from the origin, the theory has a runaway potential for the mesons 
associated with massless quarks (see, for instance, the reviews p7| , p8|). The other Nf — NfQ 
flavors are massive but much lighter than the strong coupling scale. 
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One can analyze this theory in the IR by using the Seiberg duality [11|, which trans- 
forms the model above to an SU(N = Nt — N c ) gauge theory and matter content of a 
gauge neutral Nf x Nf meson matrix <3?^ and Nf flavors of (anti-)fundamental dual quarks 
q % (qj). The superpotential is 

W = hTr'iq^lqj) - hfi 2 Tr(^ n + $22) + non-perturbative, (2.17) 

where Tr' is taken over the N color indices, Tr is over flavor indices, and we parameterize 

/ {&u)nxN 



\ 



$12 X 

$= $ 21 ($22)(N c -N f0 )x(N c -N f0 ) Y 

\ X Y (Z)N fo xN fo J 

^ 75 (x+ + X-)nxN * 

7f (P+ + P-)(N c -N f0 )xN 
\ (o")at /0 xAt J 

^ = 75^+ ~ P-)(Nc-N f0 )xN ■ ( 2 - 18 ) 
\ (°)n jo xN ) 

Note that the model in the previous subsection amounts to the case NfQ = N c — NfQ = 
N f -N c = 1. 

Again, rank conditions force us to expand around a SUSY breaking vacuum, as in (|2.4|). 
Indeed, considering the F-terms for the rank from the cubic superpotential coupling 



is at most N 



N 



f 



N c while the rank from the linear terms in the superpotential is 



Nf — NfQ. As long as NfQ < N c we cannot balance these terms and SUSY is classically 
broken. Interestingly, this condition is also necessary and sufficient for runaway behavior, 
which is induced by non-perturbative dynamics. 

So, the system settles into a SUSY breaking solution of the equations of motion, 





V 



\ 



/ fJ^NxN 



V 



(2.19) 



Expanding around this solution we discover, not surprisingly, a plethora of massive and 
massless modes, very similar to the toy model analyzed in the previous subsection. A 
notable field is, of course, the Z matrix which remains massless even after a one-loop cal- 
culation for the same reasons as in our simplified model. All the other modes are either 



-9- 



massive at tree-level or gain some positive mass squared at one-loop. An unimportant tech- 
nical difference from the toy model is that now Sy%_,5Rx_ are eaten by the supersymmetric 
Higgs mechanism. 2 

Thus, again, we need to understand the dynamics of the Z field near the origin. The 
global symmetry is SU(Nf - Nf ) x U(Nf Q ) a x U(Nf )$ x U(l) xp x U(1)r and is sponta- 
neously broken in the state ( ^19|) to SU(N C -N f0 ) x U{N f0 ) a x U(N f0 )a x U(l)' x U(l) R . 
The gauge symmetry SU(Nf — N c ) is completely Higgsed. We can use the subgroup 
U{NfQ) a x U{Nfo)v to diagonalize Z and make the eignevalues real. This simplifies the 
mass matrix along pseudo- moduli space considerably. It is actually just several copies of 
the one we considered in the previous subsection. In particular, these symmetry consider- 
ations imply that the quadratic term takes the form ~ Tr(Z^Z), so to determine its 
coefficient it is enough to turn on a single eigenvalue which is what we do in the following. 

Let us first consider the non-gauge interactions. For the purpose of the two-loop 
computation, it is straightforward to see that this model breaks up into N(N C — Nfo) 
copies of the basic interactions we considered in the simplified model. A straightforward 
way to see that is to reconsider any of the diagrams depicted in Fig. [I], e.g. the second 
diagram. There are (N c — Njq) possible Y mesons (since only one eigenvalue of Z is turned 
on) and the color of the squarks has to be matched and summed over, so we get another 
factor of N. Similar counting applies to the other two diagrams. 

Now we have to turn on gauge interactions. The basic observation here is that the 
spectrum of vector multiplets is supersymmetric over the whole moduli space ||. The 
reason is that gauge symmetry is broken in the sector of x which is decoupled in the mass 
matrix from SUSY breaking. Thus, our criteria that there has to be a p± scalar and either 
A or B particles are still applicable. In these circumstances, since gauge interactions are 
flavor diagonal, there are no vertices containing, for instance, A^,p,a. Hence, there are 
no diagrams with particles from vector multiplets which contribute to powers of Z in the 
effective potential. One could worry about new interactions between quarks from D-terms, 

2 2 

Vp = —y~] [Trq^T A q - Tr gT j4 g t ) , 

A 

where the trace is over flavor indices. In the mass basis one can see that all the interactions 
RfR^, where R\ and i?2 are real scalars, cancel. So, there are no relevant contributions 

2 However, the trace part remains massless as long as baryon symmetry is ungauged. The real part of 
the trace becomes massive via one-loop effects and the imaginary part of the trace is an exact Goldstone 
boson. 
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either from gauge interactions or from D-terms (Intuitively, we do not expect non-trivial 
effects from D-terms in the absence of accompanying fermionic loops.). 

We conclude that this gauge theory exhibits instability near the origin, with no nearby 
minimum, plausibly sloping to the runaway at large values of the Z meson. The effective 
potential in the Z direction takes the explicit form 

V = hV Ni ^~ 2 " f0) ("I " y + In 4) Tr&Z) + <D{(tfzf). (2.20) 
3. Brane Embedding 



We now embed the gauge theory of subsection 2.2 on intersecting branes in the type II A 



string theory. In subsection |3.1| we present the brane construction and review the mapping 



of its parameters to gauge theory. In subsection 3.2 we describe the perturbative brane 
dynamics - the classical forces between the branes - and its interplay with the perturbative 
quantum dynamics found in gauge theory. 

3.1 Brane Configuration 

To construct the brane configurations in type IIA it is convenient to decompose the 9 + 1 
dimensional spacetime as follows: 



d9,1 _ 10,3,1 



(3.1) 



The M 3,1 is in the directions (x° , x 1 , x 2 , x 3 ), common to all the branes. The complex planes 
C v , C w and the real line M y correspond to 

v = x 4 + ix 5 , w = x 8 + ix 9 , y = x 6 . (3-2) 

We begin with the brane configuration of Fig. §(a), whose low energy limit is the magnetic 
theory described in the previous section with \i = |l6| (for a review, see fHf). 

Fig. ||(a) presents a two dimensional slice (x,y), where x is a certain direction in v. 
The line at the bottom of the figure stands for an NS5 brane, which is stretched in the 
direction v and located at y = x 7 = w = 0. We shall call it the NS brane. The bullet stands 
for another NS5 brane, which is stretched in the direction w and located at v = x 7 = and 
y = y\ > 0. We call it the NS' brane. The x denotes a stack of Nf D6 branes, which are 
extended in the (x 7 ,w) space and located at v = and y = y 2 ; note that yi > y\. These 
are all the extended branes involved in our configurations. 
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N=N-N r D4 
r 



NS 



Figure 2: (a) is the brane construction of the magnetic theory with massless quarks, (b) describes 
its deformation by a non-zero \i parameter. In the electric langauge, A/yo flavors are still massless 
after this deformation. 



We also have D4 branes which are stretched between extended branes. There are 
N = Nf — N c D4 branes stretched between the NS and NS' branes, and Nf D4 branes 
are stretched between the NS' and D6 branes. Arrows on the D4 branes indicate their 
orientation. 

The low energy theory on the N D4 branes stretched between the fivebranes is 3+1 
dimensional M = 1 SYM with gauge group U(N). Strings stretched between these N 
"color D4 branes" and the Nf "flavor D4 branes" correspond to Nf fundamental chiral 
superfields Strings whose both ends lie on the flavor D4 branes give rise to gauge 

singlet superfields <3?j . These are coupled via the superpotential 

WVag = hq%qj. (3.3) 

This magnetic theory is the Seiberg dual of U(N C ) SQCD with Nf massless flavors [ jOj] . 

The mapping between the parameters of the brane construction and the gauge theory 
is the following. The classical U(N) gauge coupling <? mag is given by 

fmag = — > ( 3 - 4 ) 

where g s and l s are the string coupling and length, respectively. The Yukawa coupling h 
is given by 

h 2 = _»!*£_. (3.5) 

V2 - yi 

Finally, the superpotential (|3.3|) has flat directions corresponding to arbitrary expectation 
values of $ while setting q = q = 0. In the brane picture, giving an expectation value 
(<K) corresponds to moving the i'th flavor D4 brane to the location Wi between the NS' 
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and D6 branes. A non-zero expectation value (<£•) gives a mass h{$l) to the quarks q l ,<ji- 
Geometrically, this corresponds to the length of a string stretched between the i'th flavor 
brane and the color branes, and hence, 

»m = ^ ■ (3.6) 

Another deformation of this brane configuration, which is the main focus of this work, 
is to displace a stack of Nf — Nfo out of the Nf D6 branes relative to the NS' brane in the 
v direction. The resulting configuration is shown in Fig. ||(b). The separation between the 
Nf — Nfo D6 branes and the NS' brane is denoted by Ax. This brane system is the one 
studied in [18, |2(| 21 1; the latter focus on the special case Nfo = 0, while here we take 



< Nfo < N c . The configuration of Fig. |2|(b) is the energetically favorable one. 

After displacing the Nf—Nfo D6 branes, only NfQ of the flavor D4 branes, the "massless 
flavor branes," stay at their original position (with respect to NS'). On the other hand, 
N = Nf — N c of the flavor branes connect to the N color branes and move with them to the 
position v = Ax, where they are stretched between the NS and the Nf — NfQ D6 branes in 
the y direction. The remaining N c — NfQ flavor D4 branes remain stretched between the 
NS' and the Nf — NfQ D6 branes, and hence are tilted in the (y, v) space. 

In the low energy gauge theory, this deformation amounts to adding to the magnetic 
theory (|3.3|) a linear superpotential giving rise to the theory studied in the previous section, 
J2.17 ). The mass parameter \i in the gauge theory is related to Ax by [20] 



2 Ax (o 7 \ 

H =— 3. (3.7) 



Note that the brane pictures are reliable if the separations between the branes are suffi- 
ciently large and the string coupling is small. We thus set g s <C 1 and y\ , y2 — yi > l s , 



but consider the physics for generic values of Ax, similar to the study in |21|. In the 
regime for which Ax > l s perturbative string theory is reliable, and we can use it to study 
some aspects of the brane dynamics. On the other hand, in the regime where Ax is too 
small, the brane pictures are misleading, since perturbative string theory is not reliable. 
In particular, for Ax <C g s l s we should use gauge dynamics at low energies. 

3.2 Brane Dynamics 

Several phenomena in gauge theory have simple analogs in the brane construction. The 
tilted branes break super symmetry. Furthermore, they can be displaced between the NS' 
and Nf — NfQ D6 branes in the w direction. This corresponds to the pseudo-moduli $22 
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in gauge theory. 3 When Ax is sufficiently large, gravitational attraction of the tilted D4 
branes to the NS brane fixes these moduli at w = 0. In gauge theory, an analogous effect 
is the stabilization of ^22 at the origin by the one-loop effective potential. Remarkably, it 



was observed |10| , 24] that in a large class of brane constructions gravitational attraction to 
the NS brane predicts phenomena which are realized in the low energy gauge theory due 
to one-loop quantum effects. 

The other pseudo-moduli in Fig. 0(b) are the Nfo deformations of the massless flavor 
branes between the NS' and the Nfo D6 branes. These correspond to the expectation 
values (Zj) in gauge theory. The location of the Nfo D4 branes in w is not fixed by an 
attraction to the NS brane, since the NS and these D4 branes are mutually BPS. Indeed, 
the one-loop effective potential in gauge theory does not fix the pseudo-moduli Z, as we 
have seen in the previous section. 

We are thus led to consider subleading effects in the perturbative string theory regime 
(The analogous effect in gauge theory is the two-loop effective potential we studied.). There 
are several effects here which play an important role. Let us first focus on the NS' brane 
in Fig. |2](b) and further concentrate on the dynamics of the end-points of the fourbranes 
ending on it. 4 These are codimension-two objects in the world-volume theory of a type IIA 
fivebrane. To understand their interactions we can consider the effective theory in three 
space-time dimensions. In this theory the end-points of D4 branes correspond to localized 



sources giving rise to an electric field and some scalar fields, as in [29|. 

More specifically, a single D4 brane ending on an NS' brane at wq and going out in 
the direction y gives rise, for large \w — Wq\, to the following fields (in the normalization 



3 More precisely, to the expectation values ((3?22)1); non-diagonal expectation values of "1> can be seen in 
the brane pictures if one separates the D6 branes in the y direction. 

4 We thank David Kutasov for pointing out the importance of the end-points dynamics, and for very 
helpful and interesting discussions. 

5 Recalling that the world- volume theory of a fivebrane in IIA string theory does not contain vector fields 
it may be confusing that it appears (sourced by the end-point) after dimensional reduction. The point is 
that the six-dimensional theory contains five scalars. Four are the usual Goldstone modes which encode 
the shape of the fivebrane in ten dimension and the fifth is a compact scalar which has to do with the 
M-theory circle. This compact scalar has a monodromy around the D4 end-point. Upon reducing to three 
dimensions we can use Poincare duality and turn this vortex source into a usual local electric source for an 
Abelian three-dimensional gauge field. For a related analysis see |j(J. We are grateful to Ofer Aharony for 
very helpful and interesting discussions. 
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of 0; 



y = g s l s \n \w - w \, A = — In \w - w \. (3.8) 



A fourbrane going out at an angle 9 (like the tilted D4 branes in Fig. §(b)) from w' has 
the following profile (as follows by rotational invariance in the xy plane): 

V = 9sh cos 9 In \w — w' \, x = g s l s sin 9 In \w — w' \, Aq = —ln\w — w' \. (3.9) 

Since the world-volume theory of a single fivebrane is free we can use the superposition 
principle to construct a solution for two such D4 branes, 

y = g s l s (ln \w — wq\ + cos#ln \w — w \), x = g s l s sin#ln \w — w' \, 

A = -(ln\w -w \+ln\w-w' \). (3.10) 

It is straightforward to compute the binding energy of the system. Of course, scalars of 
like charges attract while identical electric charges repel. When 9 = the system is BPS 
and the forces conspire to cancel. A non-zero relative angle does not affect the electrostatic 



force, as is evident in ( 3.10 ), but it decreases the attractive force from the exchange of y 
bosons. There is no overlap in the x direction so there is no binding force from exchanges 
of x. Hence, the end-points repel. The magnitude of the repelling force behaves, for large 
\Aw\, like 

F(Aw) = "• 2( ;-. C ^ ) , (3.11) 
L s \Aw\ 

where Aw = wo — w' . 

Evidently, there are other forces acting in the system. In general, separated non- 
parallel D branes in flat space always attract since gravity dominates the RR repulsion. 
So, far away from the NS' brane our Nfo D4 branes may feel some attraction. However, 
the dominant effect near the NS' brane is expected to be the Coulomb repulsion ( 3,ll| ). 



To understand better the dynamics of this system one should solve the full non-linear 
DBI action (which should shed light on the short distance modifications of this Coulomb 
repulsion) as well as analyzing better the closed string interactions involved. Nevertheless, 
the considerations above strongly suggest that the end-points repel each other and the 
origin at w = is destabilized. 

To recapitulate, we presented an argument that the Njq D4 branes in Fig. §(b) are 
destabilized in the brane dynamics regime, nicely matching the field theory expectations. 
The analysis in the perturbative brane regime is straightforward and transparent compared 
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to the intricate two-loop computation needed in the gauge theory. However, the classical 
analysis above is not complete, but only presents some evidence for what appears to be the 
correct dynamics. It will be nice to perform more complete analysis of the various effects 
we described above and to obtain quantitative predictions for the fate of this system for 
generic separations of the brane. 

4. Comments on General Distributions of Masses 

There are some detailed implications of the results we obtained in the previous sections, 
but here we restrict ourselves to some qualitative features and postpone the complete 
phenomenological analysis to the future. Consider massive SQCD with Nf quarks in the 
free magnetic phase ordered as 

< mi < mi... < mjVp (4-1) 

where we take the mass matrix to be diagonal with positive real eigenvalues m». We are 
interested in estimating how large should the hierarchy be, and among which masses, such 
that the model is destabilized. 



By the Seiberg duality we arrive at the theory of subsection 2.2 with the superpotential 

W = hq$q-hY / f4$ ii , (4.2) 

i=i 

where fif = m^A and A is a strong coupling scale. If [i\ = then <3?n is not lifted at 
one- loop, as we have seen. Thus, the one- loop mass of $n must be proportional to On 
the other hand, there is a non-vanishing two-loop contribution. We know that it must be 
proportional to a combination of fj,\, .., [j>n c for the simple reason that if they all vanish the 
minimum is supersymmetric and the two-loop contribution vanishes. The most dominant 
two-loop contribution comes from /U/v c . 

We conclude that what is expected to affect the question of stability is primarily the 
ratio of fij\f c and The suppressing factor is, naively, the inverse loop expansion 
parameter. The correct suppression factor is supposedly even smaller due to the loop 
coefficients we calculated. 6 

In the perturbative brane dynamics regime a similar conclusion is made by very geo- 
metric and explicit means. The theory with general masses in the magnetic description has 



6 For non-zero H\ we expect the two-loop result to contain dependence on the renormalization scale which 
can render a precise estimate more complicated. 
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N c N 




NS' 



NS 

Figure 3: The embedding of ISS with general masses into string theory. The N — Nf — N c heavier 
flavors correspond to the vertical D4 branes, while the iV c lighter flavors correspond to the tilted 
D4 branes, one of whose ends lie on the NS' brane. As the mass of a light flavor is decreased, the 
D4 approaches a vertical line. 



a brane embedding shown in Fig. ||. Now there are two competing forces: the attraction 
of the tilted branes to the NS and the repulsion among them. We shall call the tilted D4 
brane corresponding to the i'th flavor the \i{ brane, % = 1, . . . ,N C . The strength of the 
gravitational attraction of the \ii brane to the NS is dictated by /Xj, as follows from eq. 
( |3.7| ). Hence, the \x\ brane experiences the smallest attraction to the NS. On the other 
hand, it is repelled from the other branes ending on the NS'. The largest repulsion is due to 
its interaction with the [in c brane, since the angle between them is the largest. It follows 



from eqs. (3.5,3.7, 3.11 ) that the strength of this repulsion is dictated by hfi^ c - Thus, 
qualitatively, we see the same behavior as in the gauge theory: the stability of the brane 
configuration is dictated by the ratio of hfj,jsr c and 

5. Summary 

We have analyzed several aspects of SU(N C ) SQCD with Nf — NfQ massive flavors and 
NfQ massless flavors in the range < NfQ < N c < Nf < 3N c /2. Through a two-loop 
computation (which was made feasible using some simplifying observations and proper 
account of remnants of supersymmetric non-renormalization theorems), we found that the 
field theory is in a runaway phase with no meta-stable states near the origin of field space. 
We have also emphasized that our results may be important for model building inspired 
by ISS like scenarios, since one is often forced to make some hierarchy of masses to take 
care of the longevity - gaugino masses tension (or to fix some other phenomenological 
problems, e.g. Landau poles). As we have shown, the meta-stable minimum in massive 
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SQCD with hierarchical masses may be destabilized due to two-loop radiative corrections. 
It will be interesting to check what constraints are imposed by solving the above mentioned 
phenomenological problems without inducing instability. 

A similar picture was obtained for the brane embedding of this model, though by 
much more elementary means. This provides an impressive test of the, yet mysterious, 
correspondence between the brane dynamics and gauge theory in SUSY breaking configu- 
rations. The brane dynamics can be applied to other systems, leading to new non-trivial 
"predictions" in gauge theory. For example, for a general mass distribution, as in section 
there are various brane predictions regarding the two-loop results in gauge theory, and it 
will be nice to test them. 
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A. Two-Loop Effective Potential 

In this appendix we briefly review some of the results of the calculation of the two-loop 
effective potential |2(|, which are used in this work. As we show in subsection |2.2| , the 
effects of gauging are irrelevant for our calculations in this paper. Hence, for simplicity we 
shall review here interacting theories of scalars and fermions. 

Consider a model with a set of real scalars Ri and Weyl fermions ipi. The masses of 
these are given by 

Cmass = --(m^RiRj - l -M lJ ^^j + c.c. . (A.l) 
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SSS FFS FFS SS 

Figure 4: Four possible diagrams which contribute to the two-loop effective potential of a theory 
with interacting scalars (represented by dashed lines) and fermions (solid lines). 



We consider a basis where the mass-squared matrices mfj and Mfj = M\ k Mkj are already 
diagonal, with eigenvalues mf and mf, respectively. Note that the (symmetric) fermionic 
matrix Mjj is not necessarily diagonal. 

The only possible renormalizable interactions in this theory are cubic and quartic 
interactions for the scalars and Yukawa interactions of two fermions and a scalar. Following 



the conventions of 26] we parameterize them as follows: 



, 1 , ;au „ „ „ 1 . iAUi „ „ „ „ I 1 



^ k R t RjR k - —X^RiRjRkRt - i^-Y IJk ^jR k + c.c.J . (A.2) 

Note that the couplings A and A' are real and symmetric under the interchange of each pair 
of indices. The Yukawa couplings l/jfc are symmetric under interchanges of spinor flavor 
indices I and J. 

In the perturbative regime one can expand the effective potential as 

Generically, the two-loop potential depends on the renormalization scale, Q. The four 
possible diagrams which can contribute to are depicted schematically in Fig. ||. We 
will further refer to these diagrams as SSS, FFS, FFS and SS respectively. Note that the 
diagram FFS appears since the masses of fermions are not necessarily diagonal. 

The contribution of each of these diagrams is parameterized by the following functions: 



m = -:(X ljk ) 2 hss(mlmlmt) (A.4) 



'si = ^(A^) 2 /sssK 2 , mj,ml 
4 2) = U^fssimlmj) (A.5) 



^ffs = d^ Jfc | 2 /FFsK,m 2 j, m 2 k ) (A.6) 



,(2) . 1 
'FFS - 2 

Vcg s = \Y IJk Y^ k J^ p ]^j,f^(mlm^ml) + c.c. (A.7) 
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The functions / can be expressed in terms of three functions I(x,y,z), J(x,y) and J(x) 
which are defined as 

J(x) =x(ln^-l) (A.8) 
J(x,y) = J(x)J(y) (A.9) 

H x , V, *) = 2 ( X ~ V ~ Z ) ln ^2 ln ^2 + \( y ~ X ~ Z ) ln ^2 111 ^2 + ( A - 10 ) 



i(z-x- y )ln^ln^ + 2xln^ + 2 y ln^ + 



where £ is defined by 



t , \ z + x-y-J?. z + y-x-R x y 

£(x,y,z) = R 21n ln ln-ln (A.ll 

\ 2z 2z z z 

z + x-y- R z + y-x-R vr 2 

Z_Ll9 ^J-/l2 

2z 2z 3 



with 



i? = y^ 2 + y 2 + - 2xy - 2xz - 2yz. (A.12) 
In terms of I and J, the functions / are given by 

fsss{x,y,z) = -I(x,y,z) (A.13) 

fss(x,y) = J(x,y) (A.14) 

fFFs(x,y,z) = J(x,y) - J(x,z) - J(y,z) + (x + y - z)I(x,y, z) (A.15) 

fpp S (x,y,z) = 2I(x,y,z) . (A.16) 

In our specific model we need only I functions with at least one argument vanishing, 
so we give them explicitly 



1(0, x, y) = (x-y){ U 2 (y/x) - \n(x/y) ln ±-/ + -(In — ) 2 - 



x — y 1 x . 9 7T 

-^(^ + 2/) + 2xln^ + 2yln^ -xln^ln^. 

In the case that two arguments vanish it simplifies further 

7(0, 0, x) = ~x(\n ^) 2 + 2x ln ^ - ^x - y x. 

However, the expression for 1(0, x, y) is still not very convenient since it contains terms 
which have no Taylor expansion around the point x = y, which appears commonly in our 
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expressions. This is a spurious singularity which cancels once all the terms are summed. 
To remove it once and for all we use Euler's identity for Dilogarithms 

Li 2 (x) + Li 2 (l - x) = - lnxln(l - x) + — , 

which gives 

1(0, x, y) = (x- y) {-U 2 {1 - y/x) - hx(x/y) In ^ + | (In ^ ) 2 ) 
5 x y x y 

-2^ X + y ) + 2xln Q2 +2yln Q2~ Xln Q2 ln Q2- 
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